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Abstract. Highly magnetised rapidly spinning neutron stars are widely consid- 
ered to be natural sites for acceleration of charged particles. Powerful acceleration 
mechanism due to unipolar induction is thought to operate in the magnetospheres 
of isolated neutron stars, bringing the particles to ultrarelativistic energies at the 
expense of the neutron star rotational energy, with inevitable emission of high en- 
ergy photons. 

The aim of this review is to present basic ingredients of modern models of magneto- 
spheric activity of rotation powered pulsars in the context of high-energy radiation 
from these objects. Several aspects of pulsar activity are addressed and related to 
spectacular results of pulsar observations with two major satellite missions of the 
past - CGRO and ROSAT. It is then argued that high sensitivity experiments of 
the future - GLAST, VERITAS and MAGIC - will be vital for a progress in our 
understanding of pulsar magnetospheric processes. 

In a conservative approach rotation powered pulsars are not expected to be the 
sources of UHE Cosmic Rays. However, several scenarios have been proposed re- 
cently to explain the UHECR events above the GZK limit with the help of accelera- 
tion processes in the immediate surrounding of newly born pulsars. Major features 
of these scenarios are reviewed along with references to contemporary models of 
magnetospheric activity. 



1 Introduction 

High energy radiation from various classes of galactic and extragalactic ob- 
jects has been observed for nearly 30 years. Large fraction of the galactic 
sources are associated with neutron stars: rotation powered pulsars (RPP), 
accretion powered pulsars (APP), cooling neutron stars, and soft gamma- 
ray repeaters (SGR). Rotation powered pulsars like Crab, Vela and Geminga 
have a long history of successful observations with baloon-born and satel- 
lite gamma-ray and X-ray experiments. Performance of old experiments had 
been, however, surpassed in terms of sensitivity, energy range, number of pos- 
itive detections, or photon statistics per object by the COMPTON Gamma- 
Ray Observatory (CGRO) and Rontgen Satellit (ROSAT). Spectacular re- 
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suits of observational campaigns of RPP with ROSAT and CGRO induced a 
new wave of interest in theoretical aspects of pulsar magnetospheric activity. 

Pair creation paradigm is a pivotal element in any model of magneto- 
spheric activity of RPP. Electron-positron pairs (e ± -pairs) are necessary since 
they are thought to be responsible for radio emission observed in radiopulsars 
and interpreted as the coherent curvature radiation of e plasma. Pairs can 
be produced in magnetospheric environments either via photon absorption 
in a dense field of soft photons (photon-photon collision) or via photon ab- 
sorption in a strong magnetic field. In either case a supply of high-energy 
(HE) photons is required in order to fulfill stringent threshold conditions for 
the pair creation. It it quite reasonable then to assume that not all of those 
HE-photons would be subject to absorption. On the contrary, many HE pho- 
tons will escape the magnetosphere without any attenuation. This argument 
leads us to expect that RPP (and all radiopulsars in particular) should be 
the sources of HE radiation. To make the production of HE photons possible, 
highly relativistic charged particles are to be injected into the magnetosphere. 
One may speculate that some of these particles will either retain their energy 
or regain it (under circumstances to be specified) upon escaping from the 
source. It is up to theoretical models of the RPP activity to show whether 
a rate of HE radiation and/or particles is interestingly high with respect to 
the sensitivity of recent and future HE detectors and telescopes. 

This review focuses on RPP as sources of HE photons, presenting the 
most important observational results as well as their interpretation in terms 
of basic processes expected in the magnetospheres of RPP. The interpreta- 
tion is offered by referring to a particular class of models of magnetospheric 
activity, known as polar gap models (or polar cap models). The name re- 
flects the association of accelerator (the gap) with a polar cap on neutron 
star (NS) surface. Contrary to polar gap models, outer gap models p3] pos- 
tulate the existence of accelerators located in regions where local corotation 
charge density reaches zero, close to the light cylinder. The e ± -pair creation 
occurs there either via one photon magnetic absorption (Crab-type outer 
gaps) or via photon- photon collisions (Vela type outer gaps). These models 
are relevant for both classical and millisecond pulsars with sufficiently high 
spin-down luminosity L s( j. The outer-gap accelerators cease to produce 
pairs once the pulsar crosses the death-line logP = 3.8 logP — 11.2 p^j . A 
modern version of the outer-gap accelerator, the so called 'thick gap solu- 
tion' is however able to accommodate pulsars of longer spin periods, like 
Geminga and B 1055-52, removing at the same time serious problems with the 
original model of |2(J. A detailed review of outer gap models in the context 
of HE radiation is available |l!| and therefore we'll concentrate on polar gap 
models. 

The review is organised in the following way: Sect. 2 defines basic quan- 
tities and introduces the assumptions used in pulsar physics. The status of 
X-ray and gamma-ray observations of RPP and essential features of the ra- 
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diation detected from several sources are presented in Sect. 3. Sect. 4 offers 
simple estimates of how effective a unipolar inductor (i.e. accelerator) can 
be when acting in the framework of a neutron star. Sect. 5 discusses verti- 
cal structure of some polar gap accelerators. Sect. 6 presents the properties 
of most important radiative processes induced by such inductors inside the 
RPP magnetosphere. With energetic arguments formulated in Sect. 4, Sect. 7 
addresses a question of whether newly born and fastly spinning RPP might 
lead to generation of ultrarelativistic charged particles responsible then for 
the UHECR events observed above the GZK limit. Sect. 8 emphasizes the 
anticipated role of high sensitivity HE/VHE missions of the near future in 
contributing to the physics of RPP. 

2 Basic parameters 

The aim of this section is to define basic quantities used throughout the 
review and to introduce their mutual relations. Several excellent monographs 
covering this subject in a detailed and sophisticated way are available e.g. 
[ [5i| with a critical discussion. 

The starting point are two major quantities measured for pulsars - period 
P, interpreted as a period of rotation of a neutron star, and P, its time 
derivative. Suppose, that a neutron star of radius R s and moment of inertia / 
rotates with angular velocity fl = 2ir/P which decreases in time (for whatever 
reason) at a rate fl = ~2ttP~ 2 P < 0. The rotational energy and its time 
derivative then read 

E mt = \ I f2 2 ~ 2 x 10 46 I i5 p- 2 erg (1) 
E Iot = / Q Q ~ -4 x 10 31 h 5 P- 15 p- 3 erg s" 1 (2) 

where P is in seconds, P-15 = P/10~ 15 and Z45 = I/10 45 g cm 2 . Instead of 
P ro t one uses the so called spin-down luminosity L s( j defined as 

L s d = —E rot . (3) 

The name luminosity is misleading since the carriers of the major part of 
E rot are not luminous for us: no one has ever managed to "see" them in a 
direct way by any type of detector (but see Sect. 7) and their nature is model- 
dependent, at least for the time being. Therefore we need a model for the 
spin-down of a neutron star. Let us assume that a magnetic dipole is attached 
to the center of a neutron star, with its moment inclined at angle a to the 
spin axis f2, and let the mean strength of the field at the stellar surface is 
B s . The magnetic dipole, rotating in a vacuum will emit energy at the rate 

imagn = ^3 B 2 S sin 2 a R 6 S Q A (4) 
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suggesting thus the following model of the neutron star spin-down: 

-^sd = imagn- (5) 

The quantity B s sin a can be inferred from P and P for a neutron star with 
known values of I and R s . For a large number of randomly oriented rotators 
the factor sin 2 a can be replaced with its averaged value of 2/3. 

Another model, where the dipolar radiation is replaced with a magneto- 
spheric wind of particles pi| , gives similar result as ([|) for an orthogonal 
rotator: 

L s d = Lwind ~ ~t Rs f2 (6) 

and therefore is independent of the angle a. Since there exists no observational 
support for P depending on sin a the standard approach is to apply the latter 
model to derive the strength of dipolar component of the magnetic field 

B\ 2 = 10 15 h 5 R 6 6 PP (7) 

where B X2 = B s /I0 12 G, and R 6 = i? s /10 6 cm. 

Assuming that B s does not change with time one can integrate (Q) to 
obtain the characteristic spin-down time scale r — a period of time elapsed 
since the pulsar was born with initial period P; 



P 
2P 



P 2 ' 



As long as Pi <C P, which is thought to be satisfied for all classical pulsars 
and most of millisecond pulsars, the last factor in (g) plays no role and thus 
t « P/2P. 

It is likely that neutron star magnetic fields contain high-order multipoles 
which may dominate the dipolar component at the surface level. Their rel- 
ative amplitudes as well as distribution remain, however, unknown. It will 
be assumed throughout the paper that the dipolar magnetic field is not dis- 
torted by rotational effects or presence of strong outflowing wind of particles 
(the latter effect has been recently invoked to decrease very high values of 
B s inferred from P and P for two SGRs ]39| ]; in consequence, their classifica- 
tion as magnetars became questionable). The field is therefore approximated 
with axisymmetric static dipole with field lines satisfying r sin -2 9 = i?dc in 
polar coordinates r and 9, with the dipole constant Rdc- The dipole con- 
stant at which rigid rotation reaches speed-of-light limit is called the light 
cylinder: R\ c — c/fi. All field lines which cross the light cylinder are then 
considered as open lines, and their footpoints on the stellar surface define 
two polar caps of radius R pc » R s ■ (i? s /i?i c ) 1 / 2 , where the latter factor is 
the sine function of the polar coordinate 9 for the outer rim of the polar cap: 
sin0 pc = (i? s /i? lc ) 1 /2. 
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3 Observational overview of high energy domain 

A posteriori evidence that high-energy activity of pulsars must somehow draw 
from their rotational energy E Iot = Ifi 2 /2 comes from a simple finding that, 
essentially, the success of detection of a particular pulsar in X-rays and/or 
gamma-rays was strongly correlated with its position in the lists of targets 
ranked by spin-down flux values L s( \/D 2 . 

The aim of this section is to review recent observational status of RPP 
in the high-energy (HE) domain, with X-rays included. The HE domain is 
hereafter arbitrarily defined as extending from a fraction of keV up to about 
30 GeV. Nevertheless, more emphasis is put on gamma-ray results. Gamma- 
ray detections are particularly precious since their interpretation is thought 
to be less ambiguous in comparison to X-ray detections. In the latter case 
(especially for very young objects) contributions from initial cooling, inter- 
nal friction, or other factors of a priori unknown magnitude may dominate 
the X-ray emission (Note: indeed, four pulsars - Vela, Geminga, B1055-52 
and B0656+14 - are classified as initial cooling candidates since their 
X-ray emission is dominated by a component which may be modelled by a 
blackbody emission from a NS surface). 

For more than 1500 pulsars known to date only about 35 positive de- 
tections in X-rays and no more than 9 detections in gamma-rays have been 
achieved. There are firm detections by CGRO of 7 pulsars (dubbed Seven 
Samurai) and another 2 cases classified as 'likely' detections. The gamma- 
ray sources were identified by virtue of flux pulsations with previously known 
P and P. Crab and Vela are the only pulsars seen by all three instruments 
of CGRO. No trace of pulsed signal in VHE range (300 GeV - 30 TeV) has 
been found so far for the gamma-ray pulsars ]8(| |64|] However, strong 
steady VHE emission is associated with 3 out of 9 gamma-ray pulsars. Two 
plerionic sources of the steady VHE radiation - The Crab Nebula and the 
plerion around B1706-44 - may serve as standard candles, with 'grade A' 
according to |pl| . A third plerion - around the Vela pulsar - was given 'grade 
B' in the same ranking. All 9 gamma-ray pulsars are strong X-ray emitters. 

The positions of these HE pulsars are shown in the P — P diagram of 
Fig.|l| along with positions of about 700 radio pulsars for which P values 
were available. A remarkable fact is that the location of X-ray sources does 
not correlate with the inferred strength of magnetic field B s ; at least not 
in a naively anticipated way that high-B objects would emit HE radiation, 
whereas low-B objects would not. In particular, 10 millisecond pulsars - about 
thirty percent of all millisecond pulsars (the objects with P ^ 0.01 s and 
P 10 -17 , i.e. with low B values: B s ^ 10 9 G) known to date - have been 
detected as X-ray sources. So far, millisecond pulsars eluded the detection in 
gamma rays and just upper limits have been available for a handful of them 
from EGRET observations §||. In the case of J0437-4715 the upper limit is 

1 /2 

interestingly tight - in disagreement with the empirical relation L 1 oc LJ d 
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Fig. 1. P — P diagram for Rotation Powered Pulsars. The pulsars detected exclu- 
sively in radio are indicated with dots; they are taken mostly from the data base of 
]7(| . Thirty five pulsars emitting X-rays are indicated with bullets. These include 
two objects recently discovered with RXTE: J0537-6910 in SNR N157B in LMC 
] is the fastest young pulsar known, spinning twice as fast as the Crab pulsar but 
with similar value of spin down luminosity; J1846-0258 in SNR Kes-75 with 
P — 0.32s, the highest P among all RPP and no radio counterpart so far. Seven 
bullets in circles indicate seven gamma-ray pulsars. Dashed lines correspond to con- 
stant values of the spin down luminosity L s d. The upper line (L s d « 4x 10 38 erg s _1 ) 
includes the Crab pulsar and J1846-0258, the lower one (L sd » 3 x 10 28 erg s" 1 ) 
includes J2144-3933 - the slowest (P = 8.5s) radio pulsar detected so far Jm)] . Dot- 
ted lines correspond to constant values of the dipolar component of the surface 
magnetic field as inferred from P and P. 
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(see Fig.^j) . Very recently, however, the likely detection of pulsed gamma-ray 
emission from J0218+4232 has been reported [ p7| . 

Spectral analysis for pulsars detected with ROSAT PSPC (0.1 keV to 
2.4 keV) shows that in most cases a power-law spectral model provides accept- 
able fits to the data |Q. Moreover, an intriguing empirical relation between in- 
ferred X-ray luminosity and spin down luminosity was found, Lx ~ 0.001 L sc j, 
confirming rotational origin of most of the X-ray activity. An interesting point 
is that the relation was obtained for all the sources regardless their temporal 
characteristics (about 50 % of all pulsars detected with ROSAT are unpulsed 
sources). Fig.|| presents these results in a somewhat different way and for 
slightly different values for Lx (compiled by the author). A complementary 
empirical relation was found for pulsed emission from 19 pulsars observed 
with ASCA (0.6 keV to 10 keV). Assuming opening angle of X-rays to be 
one steradian, the inferred pulsed X-ray luminosity correlates with spin-down 
luminosity as 

according to ||| . 

A similar power-law empirical relation holds for gamma-rays (cf. Fig.^j) , 
but with a different power-law index (e.g. ]78[]): 

L, ~ 10 35 (—^—) V2 erg s~\ (10) 
\10 38 erg s -1 / 

Important conclusion from Figs.|| and |^ is that neither Lx nor L 1 be- 
comes a sizable fraction of L s &. The most efficient conversion of spin-down 
luminosity into high-energy radiation is taking place for B1055-52 - the oldest 
pulsar among Seven Samurai - with L 1 ~ 0.1 L s d- Since in both, polar gap 
and outer gap models most of the energy gained by charged particles in the 
gaps is transferred to HE photons, it means that the wind of particles leaving 
the magnetosphere at the light cylinder is also unimportant energetically, i.e. 
-^wind "C i s d- Therefore, the lion's share of the spin down luminosity is prob- 
ably carried away in a form of magnetic dipole radiation L ma g n ; in terms of 
the so called magnetization parameter er it means cr > 1 at the light cylinder 
(cf. Sect.7). 

Broadband energy spectra per logarithmic energy bandwidth extending 
from radio, optical and UV, to X-rays and gamma-rays, constructed for pulsed 
phased-averaged components of Seven Samurai are particularly impressive 
[ fr8| and instructive. FigJ^ reveals substantial spectral differences among the 
objects, which became a subject of theoretical debates and speculations. 

Short spectral characteristics of all 9 gamma-ray pulsars (including 2 
'likely' sources) are given below after Jr§ for EGRET, |6§ for COMPTEL 
and p9] | for OSSE (and references therein). These instruments operating on 
board CGRO p3[ covered the following parts of HE domain: the Oriented 
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Fig. 2. X-ray luminosity versus spin-down age r for 29 out of 35 pulsars detected 
with ROSAT, ASCA and RXTE. Isotropic emission into 4-7T steradians was assumed 
for pulsed and unpulsed sources to infer Lx. The empirical relation Lx ~ 0.001 L s d 
found by Q can be rewritten as Lx oc B s ~ 2 r~ 2 . If all classical and millisecond 
pulsars were to have the surface magnetic field B B of a fixed value 10 12 G and 
3 x 10 8 G, respectively, they would follow the two dashed lines labelled with B B . 
Four filled circles are the initial cooling candidates; in increasing r these are: B0833- 
45 (Vela) Q, B0656+14 J0633+17 (Geminga) @, and B1055-52 @. In three 
cases the fitting of the data with either blackbody or power law spectral model 
was equally justified, but inferred X-ray luminosities are strongly model-dependent 
(circles connected with vertical bars); in increasing r these objects are: B0833-45 
||], B2334+61, and B0114+58 H}. 



Scintillation Spectrometer Experiment (OSSE) was operating in the energy 
range 50 keV - 10 MeV, the Imaging Compton Telescope (COMPTEL) - 
in the energy range 0.75 MeV — 30 MeV, and the Energetic Gamma Ray 
Experiment Telescope (EGRET) in the energy range 50 MeV — 30 GeV. 

1) B0531+21 (Crab) - Detected by EGRET, COMPTEL, and OSSE. Its 
gamma-ray flux consists of pulsed and unpulsed components, the latter one 
coming from Crab Nebula. The overall phase-averaged photon spectrum in 
the range between 50 keV and 10 GeV is described satisfactorily by a broken 
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Fig. 3. Gamma-ray luminosity versus spin-down luminosity for seven pulsars (filled 
dots) detected with the CGRO instruments. Opening angle of one steradian was 
assumed for the gamma-ray emission. Open triangles are the EGRET upper limits 
after for 350 objects, including seven millisecond pulsars. The filled triangle 
indicates the upper limit for J0437-4715 Note, that most of the upper limits 
are well above the maximum possible value for L 1 set by L 1 = L s d (dotted line). 
The dashed line marks the empirical relation derived for the CGRO pulsars: L 7 oc 
■£<sd ■ Solid lines show evolutionary tracks for a classical pulsar with B E — 10 12 G 
(upper line) and a millisecond pulsar with B B = 10 9 G (lower line) according to the 
phenomenological model of |6l| ], 

power-law shape diV 7 /de oc e~ a with a break at £b r — 120 keV, and photon 
power-law index a = 1.71 for e < £br, and a = 2.21 for e > £br- The energy 
flux is / 7 w 7.3 x 10~ 9 erg s _1 cm -2 . 

2) B1509-58 - Detected by COMPTEL and OSSE. The initial COMP- 
TEL detection was of marginal significance (~ 3CT-detection) in a narrow 
(0.75 — 1 MeV) energy band. However, recent analysis shows that the spec- 
trum extends to higher energies with a cutoff around 10 MeV jlfl] (Note: this 
new finding is not marked in Fig.|J). The energy flux at 1 MeV may be as high 
as 1.4 x 10~ 9 erg s _1 cm" 2 , but the COMPTEL point stands above the corre- 
sponding OSSE point by a factor of 4. The OSSE spectral fit between 50 keV 
and ~ 5 MeV with d7V 7 /de oc £~ 168 yields / 7 » 5.6 x 10" 10 erg s^cm- 2 . 
EGRET put strong upper limits for photon flux above 100 MeV and 1 GeV , 
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Fig. 4. Phase-averaged spectra for seven gamma-ray pulsars (Courtesy 
D.J. Thompson). 
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which clearly fall below simple power-law extrapolation of the OSSE spectral 
fit. That indicates a presence of spectral roll-over at several MeV, in agree- 
ment with the cutoff claimed by |l6). The pulsar has the highest inferred 
magnetic field (B s ~ 1.5 x 10 13 G ) among seven gamma-ray pulsars, an es- 
sential point for explaining the cutoff at 10 MeV as due to photon-splitting 
effect Q. 

3) B0833-45 (Vela) - Detected by EGRET, COMPTEL, and OSSE. Its 
phase-averaged photon spectrum between 30 MeV and 2 GeV can be re- 
produced as a power law with a — 1.7, and a strong spectral break above 
~ 4 GeV. The spectrum flattens out in the OSSE range with a = 1.3. Esti- 
mated energy flux reaches ~ 9 x 10 -9 erg s _1 cm -2 (the brightest object in 
the gamma-ray sky). 

4) B1706-44 - Young Vela-like pulsar, detected by EGRET. The spectrum 
extends from 50 MeV beyond 10 GeV and may be approximated with a 
broken power law, with photon index a changing from 1.27 to 2.25 at 1 GeV. 

5) B1951+32 - Detected by EGRET and COMPTEL. The spectrum ex- 
tends from 0.75 MeV up to 30 GeV and may be approximated with a single 
power law with photon index a = 1.89. Extremely sharp cut-off, with no 
apparent decline in the flux level. But extrapolation towards TeV falls 2 or- 
der of magnitude above an upper limit (not shown in Fig.|J) set by Whipple 
group. 

6) J0633+17 (Geminga) - Confirmed detection by EGRET only. The 
photon spectrum may be approximated with a single power law, with a = 
1.50, extending from 30 MeV to a roll-off at 2 GeV. 

7) B1055-52 - Detected by EGRET above 70 MeV. Its spectrum can be 
represented by a single power law with photon index a — 1.73, and a possible 
break around 1 GeV. 

8) B0656+14 - A 3er-detection by EGRET had been reported. The pulsar 
with its parameters (P, P, and particularly L s d) resembles Geminga nad 
B1055-52. Its photon spectrum estimated for low number of events may be 
represented between 10 MeVand 10 GeV as a very steep power-law with the 
index a = 2.8. 

9) J0218+4232 - Marginal detection by EGRET (at 3.5 a level) of pulsed 
emission has been reported recently in the energy range 100 MeV — 300 MeV 
for this distant (D <; 5.85 kpc) millisecond (P = 0.0023 s) pulsar in a bi- 
nary system with a low-mass white dwarf |47|| . The inferred luminosity of 
the pulsed emission for 1 steradian opening angle reaches L 1 ~ 1.64 x 
10 34 erg s- 1 ~ 0.07 L sd . 

Gamma-ray light curves differ significantly from those in X-rays, optical, 
and radio. Their most striking feature are relatively long duty cycles as well 
as phase shifts in comparison to the radio pulses. Only for the Crab pul- 
sar the peaks in gamma-rays as well as in radio wavelengths occur at the 
same rotational phases. The light-curve shapes fall into two categories. The 
Crab pulsar, Vela and Geminga show two sharp pulses separated in phase by 
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0.4— 0.5 and connected by an interpulse bridge of considerable level. B1706- 
44 shows two peaks separated by 0.2 in phase, with some hints of a third 
component in between. Other pulsars exhibit broad single pulses. Unknown 
opening angles for gamma-ray emission introduce a factor of uncertainty when 
inferring the gamma-ray luminosities. Broad peaks in gamma-ray pulses do 
not necessarily mean large opening angles for gamma-ray emission. Polar cap 
models, which rely on purely dipolar magnetic fields postulate nearly aligned 
rotators, where inclination of magnetic axis to spin axis is comparable to the 



angular extent of the polar cap 1 27 



With a dicovery of sharply peaked pulsed X-ray emission in the fastest 
millisecond pulsar B1937+21 an apparently separate group of millisecond 
X-ray pulsars emerges, with its members - B1821-24 @, J0218+4232 @, 
and B1937+21 |7j| - being scaled-down versions of the Crab pulsar as far 
as sharp pulse profiles and hard power-law X-ray spectra are concerned. An 
astonishing common feature within the group is the same strength of the 
magnetic field estimated at the light cylinder and the fact that it matches 
the strength of the Crab pulsar magnetic field at the light cylinder. 



4 Unipolar induction — a toy model 
4.1 Vacuum rotator 

Let us begin with a frequently invoked order-of-magnitude estimate adver- 
tising rotating neutron stars as potentially powerful accelerators and thus 
good candidates to explain UHECR (the problem addressed in Sect. 7). Con- 
sider a neutron star of radius R s and surface magnetic field B s as a perfect 
conductor. For the rotating star with its dipolar magnetic field immersed in 
a vacuum an external quadrupole electric field £ developes, with non-zero 
component along magnetic field lines at the surface. The corresponding elec- 
trostatic potential <P in polar coordinates r and 9 reads 

<P{r,9) = (3cos 2 <9- l) , (11) 

where Q tx B s /P is the quadrupole moment. Maximal electromotive force 
will then be induced between one of the two poles of the star and its equator: 

^equator - ^ B(2Rl (12) 
The corresponding voltage drop 

^equator « 3 X 10 16 B 12 R^P- 1 Volt (13) 

reaches huge values. If such a unipolar inductor was to operate in the Crab 
pulsar, B1509 or J1846-0258 (see Fig.Q), it would bring a fully ionized atom of 
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iron (Z = 26) close to the energy of 10 20 eV. However, the assumption about 
a vacuum surrounding the entire star is not correct. The space containing 
field lines closed within the light cylinder is expected to fill in quickly with 
trapped charged particles supporting the electric field which forces them to 
corotate with the star (cf. the subsection below). Therefore, the only regions 
on the stellar surface appropriate for the unipolar induction to act are those 
containing open field lines only, i.e. the polar caps. From ([II]) and ( |l2| ) the 
potential difference between the pole and the outer rim of the polar cap 
follows as 

A$ pc = Equator (^j , (14) 

with voltage drop of 

AV pc « 7 x 10 12 B 12 P- 2 Volt. (15) 



4.2 Rotating magnetosphere 

We will concentrate hereafter on a class of models where the supply of charged 
particles from a neutron star surface along open filed lines is not limited by 
binding or cohesive energy of the particles and therefore can reach the so- 
called Goldreich- Julian rate at the surface. Such a supply of charges was 
dubbed 'Space Charge Limited Flow' JFH , Q or 'free emission'. A concise 
but to-the-point account of essential properties of the SCLF models has been 
presented recently by || . 

Three boundary conditions essential for the electrodynamics above the 
polar cap are 

1) £ ■ B — for the magnetosphere within the closed field lines, 

2) (p = at the surface and at the interface between the closed magneto- 
sphere and the open field lines, 

3) £ II = at the surface level. 

Last but not least, it is assumed that the outflow is stationary and the mag- 
netosphere remains axisymmetric. 

The electric field £ required to bring a charged particle into corotation is 



-((!? - w LT ) x r) x B = 



(16) 



where the inertial- frame dragging effect is included |5^| with wlt = Kg(Rs/r) 3 J7, 
and K g = 2GI/Rlc 2 . 

Charge density necessary to support this local £ is 



J_ fl B 

Qcorot — . V • C ~ 

47T Z7TC 



1 — Ka 



(17) 
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The charge density g CO rot due to SCLF at r = R s is called Goldreich-Julian 
charge density and it is labelled with 'GJ': 



£>GJ 



n b 

2-kc 



[1 



(18) 



As charged particles flow out along the open field lines a deviation of the 
local charge density gi OC ai from the local corotation density g CO rot develops. 
By using now two relations satisfied in the dipolar structure, B(r) oc r~ 3 
and g(r) oc r~ 3 , one obtains a simple formula for local deviation from the 
corotation charge density: 



n b 

f?local Qcorot ~ 

ZTTC 



1 



Accordingly, the accelerating potential drop in SCLF reads 

3" 



A<P pc K g 



r 



(19) 



(20) 



This is a remarkable result obtained by |M : due to the inertial-frame dragging 
effect the particles drop through the potential which is significantly larger 
than in the classical approach of [0 (for P — 1 s it becomes 10 times larger). 
Moreover, the electric field £u developes now along all open lines regardless 
their orientation with respect to the spin axis (the effect is not shown in 
this simplified presentation). Therefore, all open field lines are 'favourable', 
in contrast to H. 



5 Electric field structure in SCLF gaps 

The model with frame dragging effects HU, presented in a simplified form 
in subsection 4.2, does not take into account possible feed-back effect due 
to e ± -pairs formed via photon absorption within open magnetic field lines. 
Copious pair formation occurs in a relatively thin layer called for this reason 
a pair formation front (PFF). The creation of pairs leads to screening of the 
accelerating field £ m within the layer of PFF. A detailed picture of this effect 
would require to follow the dynamics of electrons and positrons in a self- 
consistent way. Instead, it is reasonable to assume, that the field is shorted 
out at the height were the first e ± -pair is created (hereafter denoted as h c ): 
£n = for h > h c . 

The problem of electric field structure in the context of SCLF with bound- 
ary condition £h = set at ho = (stellar surface) and at h > h c (PFF) was 
formulated and solved by pCfl . The solution is rather lengthy and includes 
special functions. It is however possible to obtain simple but quite accurate 
analytical approximations. As long as the length h c of the accelerator is of 
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the order of polar cap radius R pc , the accelerating electric field may be ap- 



where B>n = £> S /10 12 G, P is the spin period in seconds, x is the angle 
between the spin axis and the magnetic moment of the rotating star, h is 
expressed in cm, and M = XAMq, R s = 10 6 cm. The magnetic colatitude 
£ = 6/0(rj) is scaled with the half-opening angle of the polar magnetic flux 
tube 9(r/), where r\ = 1 + h/R s . The magnetic colatitude function is a 

monotonically decreasing function, with /i(0) ~ 1 and /i(l) = 0. 

Vertical structure of the electric field depends (via the location of PFF) 
on radiative processes which induce the pair creation: curvature radiation 
and inverse Compton scattering on soft X-ray photons from the stellar sur- 
face (brief characteristics of these processes is presented in the next section) . 
An interesting effect was noticed in this context |l0|]: Suppose that a small 
fraction of positrons is stopped by a residual (non-zero) electric field at the 
site of their creation and then forced to flow towards the stellar surface (the 
effect noticed already by |7^]). The backflowing positrons are expected to 
induce a formation of an additional PFF, which would short out the electric 
field at ho- These positrons cool upon the action of ICS and CR. With rea- 
sonable surface temperatures (T s ~ 5 x 10 5 K) it is ICS which dominates the 
cooling of upward moving electrons and downward moving positrons. There- 
fore, the e ± -pair creation will be induced by upscattered photons from the 
stellar surface, rather than by curvature photons. The situation is not sym- 
metrical, however, for electrons and positrons since the field of soft photons 
is not symmetrical with respect to both types of particles. In consequence, 
the positrons cool more efficiently, the ICS-induced cascades are easier to 
achieve for positrons than electrons and thus a lower PFF (where E\\ = 0) 
tends to be located above the stellar surface. Such a situation is not stable, 
therefore. However, elevating the accelerator up to altitude ho ~ 1 Rs above 
the surface diminishes the role of the ICS; the CR cooling dominates here 
and a stable accelerator is possible. No self-consistent calculations of such a 
'sandwich-like' accelerator exist at present, but the results obtained by jl0| 
with an approximate treatment of the problem look promising indeed. 

6 Radiative processes in pulsar magnetospheres 

Cooling of ultrarelativistic electrons via curvature radiation (CR) and mag- 
netic inverse Compton scattering (ICS) are the most natural ways of produc- 
ing hard gamma-rays capable of inducing cascades of e ± -pairs and secondary 
HE photons. These two processes dominate within two distinct ranges of 
Lorentz factors 7 of primary electrons. 



proximated according to as 




(21) 
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When 7 ^ 10 6 , magnetic inverse Compton scattering plays a dominant 
role in braking the electrons and it is the main source of hard gamma-ray 



photons 1 73 . Energy losses due to resonant ICS limit the Lorentz factors of 
the particles to a level which depends on electric field strength Sj, tempera- 
ture T and size of hot polar cap, and magnetic field strength B s JF^]. The 
Lorentz factors can then be limited even to ~ 10 3 . This stopping effect be- 
comes more efficient for stronger magnetic fields, and it was suggested as an 
explanation for the observed cutoff at ~ 10 MeV in the spectrum of B1509-58 
& 

However, in their modern versions the accelerators of particles are strong 
enough to outpower the ICS cooling. In consequence, very high Lorentz fac- 
tors - 7 10 6 - are achieved by electrons, limited by CR. The first detailed 
scenario of radiative processes in CR-induced cascades was presented by Q 
and despite many modifications and additions its basic features remain valid. 
The model assumes that primary electrons accelerated to ultrarclativistic en- 
ergies emit curvature photons which in turn are absorbed by the magnetic 
field and e ± -pairs are created. These pairs cool off instantly via synchrotron 
radiation process. Whenever the SR photons are energetic enough they may 
lead to further creation of pairs, etc.. ICS can still be incorporated to the 
models with CR-induced cascades as the process relevant for e ± -pairs, since 
typical Lorentz factors of theirs do not exceed ~ 10 3 . According to the an- 
alytical model of |85| the empirical relations for X-ray and gamma-ray lu- 
minosities of pulsars (presented in Sect. 3) can be reproduced satisfactorily 
when the ICS involving e ± -pairs is included. 

Processes relevant for production and transfer of HE radiation in pulsar 
magnetospheres are, therefore: 

- Curvature Radiation, 

- Magnetic Inverse Compton Scattering, 

- Magnetic Pair Creation 7 — > e 

- Synchrotron Radiation, 

- Photon splitting 7 — > 7 + 7, 

- Photon-photon Pair Creation 7 + 7 — > e ± . 

Basic properties of these processes are briefly reviewed below. The last 
process in the list has been omitted. The reason is that, whenever recalled in 
the context of pulsar magnetospheres, photon-photon pair creation is treated 
exactly as in free space. Such treatment is justified in models of 'thick outer 
gaps' but within the framework of polar cap models this is not the case, 
in general. However, no handy formula is available for the cross section of 
this process in the limit of high B and standard non-magnetic formulae are 
in use whenever necessary. 

In order to illustrate the significance of these processes in forming HE 
spectrum extending over many decades in energy, the numerically calculated 
effects due to the first 4 processes in the list will be presented in Fig.|| (after 
[p0| ) along with overlaid data points for the Vela pulsar. The dipolar field 
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in the Vela pulsar does not exceed 10 13 G and thus photon splitting is not 
competitive to magnetic pair creation; its effects are negligible and will not 
shown. The electric field structure of the accelerator used in these calculations 
is taken after Q and appropriate rescaling due to ho > 0. However, the case 
calculated for Fig.|| was chosen with ho = 3R S , i.e. much higher than in [ ftp] ], 
in order to better reproduce the phase averaged spectrum of the Vela pulsar. 



6.1 Curvature radiation 

Relativistic electron of energy 7771c 2 (we take 7>1) sliding along the mag- 
netic field line of curvature g CI will emit photons with a continuum energy 
spectrum peaked at 

e poak w0.29e cr , (22) 

where 

3 7 3 

£« = rcS-, (23) 

2 Qcr 

is called the characteristic energy of CR. The radius g CI for a purely dipolar 
line attached to the outer rim of the polar cap can be approximated not far 
away from the NS surface as g„ ~ %/ Re ■ Ric ~ 10 s s/Pcm. The cooling rate 
of that electron is 

2 e 2 -y 4 

7cr = -~ — \- (24) 

For a monoenergetic injection function of electrons (5(7) oc £(7 — 70) and 
their cooling due solely to CR the electrons will assume a single power-law 
distribution in energy space AL(el.) oc 7 -4 for 7 < 70 as long as they stay 
within the region of the cooling. Their escape introduces a natural low-energy 
cutoff 7cutoff in Nj(el.). Therefore, the unabsorbed CR energy spectrum f e (e) 
due to the injected electrons has a broken power-law shape, with a high- 
energy limit set by 70 and the break at some energy Ebreak- For e > ffbreak 
the energy spectrum is / e (e) oc e~ 2 / 3 , and / e (e) oc e +1 / 3 for e < ebreak- 
Since nonthermal spectra cover ususally many decades in energy it is more 
convenient to use ef e (e) for easy comparison of power in different parts of 
energy space (see Figs|| and ||). Accordingly, ef £ (e) oc e +1 / 3 above the break, 
and oc e +4 / 3 below the break. 

The cutoff limit 7 cu toff can be found by comparing the characteristic cool- 
ing time scale t CT = 7/|7 cr | with the estimated time of escape 7j csc , which we 
take as t csc ~ g C r/c. Therefore 

ebreak » 0.29 ■ ~h — « 44 MeV, (25) 

4 r 

where r is the classical electron radius. Note, that the photon energy Ebreak 
at which the spectral break occurs does not depend on any pulsar parameters. 

The spectrum of CR calculated numerically to model the Vela pulsar |30| 
is shown in Fig|| as dot-dashed line. High-energy cutoff due to one-photon 
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magnetic absorption occurs around 10 GeV. Note the importance of gamma- 
ray detectors capable to operate above 10 GeV for (in)validating the model. 
The low-energy CR spectral break £break is prominent at ~ 40 MeV. Below 
£break the power of CR decreases and eventually becomes unimportant at 
~ 1 MeV where the synchrotron component takes over. 



6.2 Magnetic pair creation 

Pair creation via magnetic photon absorption (7 + B — > + B) is kine- 
matically correct since the magnetic field can absorb momentum. To ensure 
high chances for the process to occur it is not enough for a photon propa- 
gating at some angle a to local B to satisfy the energy threshold condition, 
sini/; • e > 2mc 2 , but high optical thickness t~ ( b within the magnetosphcrc 
is required. In fact the condition r 7 s = 1 has been used as a criterium 
for the so called death-line for radiopulsars in the P ~ P diagram. Maxi- 
mal values of sin 1 ?/; for curvature photons in the dipolar field do not exceed 
~ 0.1 sin6» pc S3 0.0014 p- 1 / 2 @, so the both conditions are difficult to meet 
for long-period rotators. 

The absorption coefficient for the process as described in [||| and used to 
calculate r 7 s reads 

v(e) = ~^T( X ) (26) 

where a is a fine structure constant, A c is a Compton wavelength, _B C rit = 
m 2 c 3 1 eh ~ 4.4 x 10 13 G, P>± is the component of the magnetic field perpendic- 
ular to the photon momentum, and x = \ j~ m <•> ^ s ^ ne Erber parameter \- 
The function T(x) is then approximated as T(x) S3 0.46 exp (— 4//3x), which 
is valid for x 0.2; for x ~ 0.2 this approximation starts to overestimate ij. 
The function / is the near-threshold correction introduced by important 
in particular in the case of classical pulsars. Electron-positron pairs created 
through the magnetic absorption radiate then via the synchrotron process 
(SR) described below. 



6.3 Synchrotron radiation 

Consider a particle of energy 7771c 2 gyrating around a local field line at a pitch 
angle ip. Let 711 denotes the Lorentz factor of the reference frame comoving 
with the center of the gyration. As long as 7|| ^> 1 it relates to the pitch 
angle ip via sin ip ss 77 1 . The energy available for synchrotron emission at 
the expense of the particle is 7j_ rnc 2 , and 7 = 7_l 7||- 
The rate of SR cooling reads 



Neutron Stars 19 



In comparison to the CR cooling it is enormous (due to much smaller curva- 
ture radius). 

Critical photon energy (analogous to (H)) reads 

£sr = 7:?i 7 2 smV'. (28) 

2 m c c 

For a monoenergetic injection function of particles (e ± -pairs in the con- 
text of this review) and their cooling due to SR the energy spectrum of SR 
spreads between a high-energy limit e sr (7o) set by 70 of the injected (created) 
particles, and a low-energy turnover e ct determined by the condition 7_i_ ~ 1: 

x 3 fi5 1 
£ ct = esr7 = 7|| ) = -% — . (29) 

The spectrum assumes a single power- law shape / E (e) oc e -1 / 2 (and ac- 
cordingly - ef s (e) oc e +1 / 2 ) above the turnover. Below e ct , the spectrum 
f e changes it slope, asymptotically reaching oc e +2 . It is built up by contri- 
butions from low-energy tails emitted by particles with 7x^1; and each 
low-energy tail is assumed to cut off at local gyrofrequency, which in the 
reference frame comoving with the center of gyration is wr = — — — . 

The spectrum of SR calculated with Monte-Carlo method to model the 
Vela pulsar is shown in Fig|| as a dashed line. The low-energy part of the 
SR spectrum at e c t seems to be essential for reproducing an interpolation 
between the RXTE and the OSSE data. 



6.4 Magnetic Inverse Compton Scattering 

Consider an electron with a Lorentz factor 7 moving along a magnetic field 
line B and a photon of energy s = emc 2 moving at angle arccos \i to the field 
line. In the reference frame comoving with the electron (primed symbols) the 
counterpart of the free-space Compton formula, due to energy-momentum 
conservation appropriate for collisions with the electron at the ground Landau 
level both in the initial and final state, reads 



4 = (l-/4 2 ) 1 jl + e'(l-///4) + 



- "l + 2ey>' s - fx') + e'V* - M') 2 ] } (30) 

where e' = £7(1 — /3fi) and symbols with no subscript and with the 

subscript 's' refer to the state before the scattering and after the scattering, 
respectively. A longitudinal momentum of the electron in the electron rest 
frame changes due to recoil from zero to (e'/i' — e' s fi' s )mc. 
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Fig. 5. The model energy spectrum calculated by |k| to reproduce the spectral 
features of the Vela pulsar (P = 89 ms, B B — 6 x 10 12 G). The accelerator is 
located at ho = 3R S above the surface (see Sect. 6). The broad-band spectrum 
consists of four components due to: curvature radiation of primary electrons (dot- 
dashed), synchrotron radiation of secondary e ± -pairs (dashed), inverse Compton 
scattering of surface X-ray photons on the e ± -pairs (thin solid) and the blackbody 
surface emission (dotted). The surface temperature T B = 1.26 x 10 6 K was assumed. 
Total spectrum is given by a thick solid line. Phase-averaged data points for Vela 
from different satellite experiments are indicated. Filled squares - RXTE ; open 
circles - OSSE @; open squares - COMPTEL Q; filled circles plus upper limit 
just above 10 GeV - EGRET Vertical axis is in log of MeV cm" 2 s" 1 units. 




The polarization-averaged relativistic magnetic cross section in the Thom- 
son regime may be approximated with a nonrelativistic formula |2q] : 



ff =^(l- M ' 2 + (l + /x' 2 ) 



91 



92 - 9i 



(31) 



where a T is the Thomson cross section, and gi and <?2 are given by 



9i{u) 



(u + l) 2 ' 



92 (u) 



(u - l) 2 + a 2 



(32) 



where u = e'/e E 



2ae B /3, e B = heB /m 2 c 3 and a is a fine-structure 



constant. The resonance condition for the scattering is therefore the cyclotron 
resonance e' — e B . The factor a represents a 'natural' broadening of the 
resonance due to finite lifetime at the excited Landau level. 

In the Klein-Nishina regime (e' > 1) the relativistic magnetic cross section 
for the |/x'| ~ 1 case becomes better approximated with the well known Klein- 
Nishina relativistic nonmagnetic total cross section a KN p5| pgfl. 
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The rate 1Z of scatterings subject by an electron moving across the field 
of soft photons, measured in the lab frame is 



M (33) 



where Q = d/id<p is the total solid angle subtended by the source of soft 
photons, /i = cos0, a is a total cross section for the process, and dn p h/ds/dfi 
is the local density of the soft photons. 

The properties of the field of soft photons are usually simplified by taking 
G?n p h/ ds I ' dfl as for the blackbody radiation. This simplification should be 
taken with care since magnetised atmospheres of neutron stars introduce 
strong anisotropy as well as spectral distortions to the outgoing radiation 
[ [58| . Effectively it means that the ICS effects obtained with this simplification 
are just upper limits to the actual effects. 

To estimate electron cooling rate 7 ICS due to the ICS the differential form 
of (|3l|) is necessary: 



da 3a T 
dJK~16n 



(1-// 2 )(1- M ' S 2 ) + 



+ ~ (1 + // 2 )(1 + a4 2 )(.9i+.92) 



(34) 



(eg. |2(|), where dQ' s — d<j)' s dfj,' s is an increment of solid angle into which 
outgoing photons with energy e' s in the electron rest frame are directed. The 
mean electron energy loss rate then reads 

(^) (1- ft.) X 

<^(^>-- f > (35) 

where e s = e' s 7(l + /3/i' s ) is the scattered photon energy in the lab frame 
(e.g.g§). 

The spectrum of magnetic ICS calculated numerically to model the Vela 
pulsar is shown in Fig.|^ with thin solid line. The blackbody soft photons 
originating at the stellar surface (dotted line) are upscattered by secondary 
e ± -pairs at the expense of their "longitudinal" energy j^mc 2 , assumed to 
remain unchanged during the burst of synchrotron emission. Without the ICS 
component due to the e ± -pairs the RXTE data for the Vela pulsar would be 
difficult to reproduce within the model. It is worth to note, that the magnetic 
ICS component due to primary electrons (not shown in Fig.^) is energetically 
insignificant comparing to the CR component because the case of strong 
accelerating field was used in this particular model. 
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6.5 Photon splitting 

Photon splitting into two photons in the presence of magnetic field B is 
a third-order QED process with no energy threshold The attenuation 
coefficient, after averaging over the polarization states, reads [Q 

We) « f iiV f^sV e * cm-, (36) 

(provided B does not exceed -B CI -it substantially) where a is a fine structure 
constant, e is the photon energy in units of mc 2 , #kB is the angle between 
photon momentum vector and the local magnetic field. The process, therefore, 
strongly depends on magnetic field strength B. 

Photon splitting has attracted substantial interest in recent years due to 
the discovery of neutron stars with supercritical magnetic fields (i.e. magne- 
tars; see Fig. 6) |50| . It has been analysed in details by and incorporated 
in a Monte Carlo code tracing the propagation of electromagnetic cascades in 
the magnetospheres of high-i? pulsars. The effect was found to explain satis- 
factorily the unusual cut-off observed in the gamma-ray spectrum of B 1509-58 
(see Sect. 3). Generally, it becomes competitive to the magnetic pair creation 
for dipolar magnetospheres with B s ^ 0.3B cr i t @. The degradation of pho- 
ton energy in the course of splitting inhibits also any development electro- 
magnetic cascades. In consequence, high-i? RPP should not emit coherent 
radio emission. Indeed, there exists a high-i? region in the P — P diagram 
(Figs. 1,6) void of radiopulsars. Even though recently discovered (during The 
Parkes Multibeam Pulsar survey) two high-i? radiopulsars JItJ are located 
above the limiting line derived by 0, the general argument for magnetars 
expected to be radio-quiet RPP remains valid flS6| . 



7 RPP as sources of UHECR generated beyond the 
light cylinder 

The hypothesis that cosmic ray events above ~ 5 x 10 19 eV, i.e. above the 
GZK cutoff, are due to charged particles accelerated by strongly magnetised 
neutron stars is being kept under consideration JlJ] ^6|Jll| apparently for 
two reasons. First, within Down-Top scenarios which rely on conventional 
physics a list of classes of objects satisfying a necessary condition to generate 
such cosmic rays is rather short according to the appealing Hillas diagram 
[ [f2| . Second, no compelling breakthrough has been achieved in studying other 
candidates in this context, though central engines of AGN or jets extending 
from Fanaroff- Riley II radio galaxies |12] take rather high positions in the 
list f|. 

In the case of neutron stars the most promising and natural reservoir of 
required energy is the rotational energy of the stars, a point reiterated on 
numerous occasions, e.g.p9|, and assumed in most models. There are three 
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fundamental features discriminating the models: 

1) the site of acceleration with respect to the neutron star, 

2) the mechanism allowing to tap the rotational energy by charged particles, 

3) the nature and origin of the charged particles subject to acceleration. 

Simple but sounding arguments supported by observed HE-radiation prop- 
erties of RPP make the first point rather clear: the process of particle acceler- 
ation should take place beyond the light cylinder. Potential advantage of this 
choice over acceleration process within the light cylinder is twofold. First, 
full potential drop across open field lines A<& pc (|14|) is available, at least in 
principle, for particles outside the light cylinder, whereas the capability of an 
accelerator inside the magnetosphere is severely constrained by copious for- 
mation of electron-positron pairs which short out the electric field £ u easily, 
i.e. A<&n < A<P pc . Second, unlike inside the magnetosphere, the acceleration 
of charged particles is not limited by any radiation losses - a point espe- 
cially important in the context of UHECR. Therefore, a particle of charge 
Ze reaches then maximal possible energy 

£ max « Ze A$ pc = 6 x 10 19 Z B 13 F ms 2 eV, (37) 

(where B13 = i? s /10 13 G) which is substantially higher than the energy E p 
which can be attained in the accelerating field of either the polar gap or outer 
gap: E p < ZeA$\\ < E maK . 

The idea of UHECR events due to an accelerator or a converter of the ro- 
tational energy into kinetic energy of particles, located beyond the light cylin- 
der has been pursued recently in a couple of diametrically different models. 
Historically the first model considered the fate of charged particles injected 
into the plerionic nebula powered by a shocked relativistic wind from a cen- 
tral pulsar |l(J [H . This model was motivated by the theoretical analysis of 
magnetohydrodynamics within the Crab Nebula |60) M]. The accelera- 
tion was proposed to occur in the electric field of the shocked wind. A simple 
structure of the electric field was derived as induced by the relativistic ra- 
dial wind crossing a toroidal magnetic field (which originates inside the light 
cylinder, at the stellar surface) with assumed radial and angular dependence 
after j(5(J and ElJ]. Charged particles entering the nebula are subject to the 
£ x B drift as well as to the WB drift (due to strong inhomogeneity in B) 
along various paths, and may gain energy at the expense of £ before exiting. 
Since the electric field is potential here, the net change in the particle energy 
does not depend on the path but solely on the points of entry and exit. In par- 
ticular, the maximal gain AE of energy is due to potential difference between 
the pole (fixed by the rotation axis of the pulsar) and the equator (cf. (3) in 
[[f8|), and actually it equals to i? max given by (|37|). In order for the particle 
to penetrate the nebula from outside it should be already highly relativistic, 
with E'init ~ 10 15 eV, presumably pre-accelerated by diffusive acceleration at 
the outer shock where the supernova remnant meets the interstellar medium. 
Otherwise the Larmor radius of the particle is to too small for the VB in- 
ward drift at the pole (the case relevant for one combination of signs of B 
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and Ze) to counter the £ x B outward drift. Though very attractive, the 
model is unable to make any strong predictions about spectral properties 
of its UHECR without reliable assumptions about spatial distribution and 
directional properties of the pre-accelerated particles. 

All other models make use of charged particles, either protons or iron 
nuclei, coming from within the light cylinder i.e. supplied by a neutron star 
itself. 

A recently proposed phenomenological model of UHECR within the Galaxy 
incorporates iron Fe 26 nuclei (Z = 26 Z2&) accelerated in a relativistic MHD 
wind flowing out of very young, rapidly rotating highly magnetised neutron 
stars The number of nuclei crossing the light cylinder per unit time is 
taken equal to the Goldreich- Julian rate estimated at the polar cap 

N(Fe 26 )^N GJ =A pc p±, (38) 
Zee 

where A pc ~ nR^R^, 1 and qqj is defined in (|i~8|). It is also assumed that 
electron-positron pairs formed within the magnetosphere do not dominate 
the flow in terms of the rest mass. This is a reasonable assumption because 
otherwise the number of created pairs per nucleus would have to ex- 
ceed mFe/2m c ~ 5 x 10 4 which is unlikely in very strong magnetic fields 
B s > 10 13 G due to photon-splitting effects |3£|. A key postulate of the model 
is that transfer of a sizeble fraction (£ ^ 1) of the spindown luminosity oc- 
curs just beyond the light cylinder into the kinetic energy flux of the ions. 
At the light cylinder, however, the spindown luminosity (i.e. the rotational 
energy loss rate) is usually thought to be dominated by the Poynting flux 
(of low- frequency electromagnetic waves). In terms of the so called magne- 
tization parameter a which is defined as the ratio of the Poynting flux to 
the the particle kinetic energy flux |}0| [|T] it means that - whatever the 
reason - a does not remain constant in the wind: a 3> 1 at the light cylin- 
der converts further out to a <C 1, i.e. the wind must depart from an ideal 
MHD flow case. This picture is motivated again by the wind models for the 
Crab Nebula (e.g. jfjlf fl45f , see also Q for calorimetric properties of the 
Crab Nebula, and |3j for an account of the problem of the coupling of RPP 
to plerionic nebulae). No consensus about likely mechanisms responsible for 
the dissipation of the Poynting flux has been reached so far, though several 
models have been proposed (e.g. jH) Q). For a single average iron nucleus 
the postulated conversion means acceleration up to 

£uhecr « L sd /N(Fc 2e ) « Ze A$ pc « 10 21 Z 26 B 13 P m 2 eV, (39) 

where (Q) was used. Similarly as in the previous model, the energy gain is in 
fact directly related to the full potential drop A<P pc across open field lines. For 
a fully ionized iron (2^6 = 1) to reach Exjhecr = 10 20 eV the pulsar would 
have to be fastly spinning and highly magnetised (for example P = 0.001s 
and B s = 10 12 G). According to p5| , the region in the P — B s space occupied 
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Fig. 6. The high-B s part of the P — P diagram (see Fig.|l]). Dashed line indicates 
the full potential drop A$ pc across open field lines equal to ^ x 10 20 V. The allowed 
region in the model of jEiJ to accelerate a fully ionized atom of iron (Z — 26) to the 
energy of 10 20 eV and let it traverse the pre-supernova envelope without spallation 
effects lies to the left of two solid lines which are drawn for two values of the mass 
of the envelope: 5M and 50Mq. A group of magnetars, clustered around P ~ 10s 
is also indicated after jio) although their X-ray activity is not driven by rotational 
energy losses (i.e. they do not belong to the class of RPP) but rather involves 
accretion or the decay of strong magnetic field. Open squares and crosses denote 
SGR and AXP, respectively. 
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by pulsars satisfying the requirement Stjhecr = 10 20 eV is actually even more 
constrained then ( p9| ) suggests. The magnetic field of the pulsar should not be 
too strong to make sure that the characteristic time scale of spin-down r is not 
too short allowing thus the expanding pre-supernova envelope of mass M cnv 
to disperse and become transparent so the iron nuclei would escape it without 
significant spallation. For M anv equal 5M Q and 50Af Q the magnetic field B s 
should not exceed ~ 6 x 10 13 G and ~6x 10 14 G, respectively (see Fig.||). A 
factor not discussed in |]l5| but likely to further constrain the allowed region in 
the B s — P space in order to make the model work is a poorly known cohesive 
or binding energy Ae c of iron on the high-_B s neutron star surface. For B s > 
10 13 G Ae c may exceed 5 keV [[[]. Thermionic emission of iron nuclei is then 
impossible unless the surface temperature exceeds ~ 3 x 10 6 K within the first 
~ 10 8 s [63 after the neutron star formation which is a characteristic time 
scale mediated by the neutron star's spin-down when the energy E max drops 
below a required CR energy ~ 10 20 eV, and the expansion of the envelope 
during which it becomes transparent to the Fe nuclei. The only supply of 
iron nuclei may then occur via field emission. Making modest assumptions 
about the distribution of rapidly rotating neutron stars, the rates of their 
birth and - most importantly - the efficiency rate for acceleration of the 
particles at the light cylinder, an expected flux of these particles on Earth 
is compared to the data. The calculated particle spectrum above 5 x 10 19 eV 
is flat, iVpo oc with 7=1, and therefore should be visible as a distinct 
component in the CR spectrum, which is much steeper below 5 x 10 19 eV with 
the power-law index 7 rj 3. However, the number of events above 5 x 10 19 eV 
detected so far is insufficient for describing them in terms of spectral features 

For a pulsar capable of producing UHECR as in [|l5| , a natural continua- 
tion - as it slows down its rotation and becomes a Crab-like object - would 
be to generate CR particles of lower energy. Hadronic model proposed by 
[|j seems to be relevant in this context, since it describes the fate of iron 
nuclei extracted from the pulsar surface and then subject to acceleration. 
However, a different acceleration site is considered here: the nuclei are as- 
sumed to accelerate within the light cylinder, in outer gaps of [^o) (i.e. the 
type of accelerator not addressed in this review). Once the nuclei achieve high 
energy (7f c > 10 5 ) they disintegrate in the field of soft photons present in 
the outer gap. Relativistic neutrons extracted in this way decay into protons 
cither inside or outside the surrounding nebula. In the first case, these rela- 
tivistic protons interact with nebular material via pp processes giving rise to 
VHE photons and neutrinos. In the second case, cosmic ray protons peaking 
at ~ 10 15 — 10 16 eV (close to the knee in the CR spectrum) are produced. 

A distinct model of UHECR particles above the GZK cutoff has been 
recently proposed J3(J . It postulates that UHECR are protons accelerated in 
reconnection sites above the magnetospheres of very young pulsars formed 
by accretion-induced collapse of white dwarves (AIC). Promising candidates 
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would occupy roughly the same region of the P — P space as indicated for the 
model of (the less massive envelope case) in Fig.[| Since the estimated 
rate of AIC events per galaxy is too low to rely just on our Galaxy, the 
contribution of all galaxies within the distance of about 50 Mpc is necessary 
to reach the observed rate of UHECR events. 




Fig. 7. Cummulative spectral flux of photons expected for the millisecond pulsar 
J0437-4715 at a distance of 140 pc according to p6| . The shaded region shows the 
range of flux levels due to uncertainity in the maximal energy of primary electrons. 
The main part of the spectrum is due to curvature radiation of the electrons. The 
additional feature reaching the VHE domain is due to inverse Compton scattering 
of soft photons from the surface with the temperature 4 x 10 K. Sensitivities of 
EGRET as well as three major HE and VHE experiments of the future are also 
indicated. MAGIC-LZA denotes sensitivity of MAGIC in its Large Zenith Angle 
mode. 



8 Concluding remarks 

High energy astrophysics of neutron stars received an impressive boost from 
two major satellite missions of the past - ROSAT and CGRO, backed by 
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still ongoing experiments - ASCA, BeppoSAX and RXTE. Theoretical as- 
trophysics of RPP was confronted with- and surprised by an unprecedented 
variety of spectral and temporal properties among the detected sources. An- 
other unexpected challenge came from radio-astronomy, due to superb per- 
formance of The Parkes Multibeam Pulsar Survey, with recent discoveries 
of radiopulsars, with extremely high magnetic fields in two cases and an 
extremely long spin period in another case. Numerous modifications (both, 
minor and major) to the existing models of magnetospheric activity were 
invented to accommodate at least some of these properties. Several predic- 
tions have been presented which would hopefully discriminate between those 
models. It will be impossible, however, to verify those predictions without 
achieving higher sensitivities and exploring new energy domains. 

A break-through in understanding rotation powered pulsars and their ple- 
rionic environments should then come from HE/VHE astronomy of the near 
future, with its planned satellite and ground-based experiments. Expected 
sensitivity and energy range for some of them is presented in Fig.^ along 
with a predicted flux from a nearby millisecond pulsar, overlaid for the sake 
of comparison. The satellite experiment GLAST [Q will be superior to its 
predecessor - EGRET on board CGRO - in two aspects. First, its sensitivity 
at 10 GeV will be three orders of magnitude better than that of EGRET. 
Second, it will reach energy of 300 GeV, closing thus for the first time a wide 
gap in energy between ground-based and satelite experiments. The MAGIC 
Telescope |L3| - a 17 m diameter Imaging Air Cherenkov Telescope (IACT) 
- is expected to operate with sensitivity about three orders of magnitude 
higher at 10 GeV than EGRET. Its advanced technology will make possible 
to cover energy range between 10 GeV and 1 TeV, and to reach ~ 50TeV 
in the Large Zenith Angle mode. Energy ranges of GLAST and MAGIC will 
overlap over more than one decade in energy. Another proposed IACT, VER- 
ITAS p^ |, will be an array of seven 10 m telescopes, covering energy range 
from 50 GeV to 50 TeV with planned sensitivity at 1 TeV about ten times 
better than MAGIC. 

Acknowledgments 

I am grateful to J. Arons, W. Bednarek, T. Bulik, K.S. Cheng, J. Dyks, 
E.M. de Gouveia Dal Pino, M. Gros, P. Haensel, O.C. de Jager, A.V. Olinto, 
S. Shibata, T.C. Weekes, and B. Zhang for sharing with me their views and 
opinions on issues addressed in this review. Financial support by KBN grant 
2P03D02117 is acknowledged. 

References 

1. A.M. Abrahams, S.L. Shapiro: ApJ 374, 652 (1991) V.V. Usov, D.B. Melrose: 
Aus.J.Phys. 48, 571 (1995) 



Neutron Stars 29 



2. S.L. Adler: Ann.Phys. 67, 599 (1971) 

3. J. Arons: Mem. Soc. Ast. Ital. 69, 989, (1998) 

4. J. Arons, E.T. Scharlemann: ApJ 231, 854(1979) 

5. J. Arons: 'Pulsar Death at an Advanced Age'. In: IAU Coll. 177, Pulsar 
Astronomy - 2000 and beyond ed. M. Kramer, N. Wex, R. Wielebinski (ASP 
Conference Series 202: 2000) pp.449-454 

6. M.G. Baring, A.K. Harding: ApJ 507, L55 (1998) 

7. W. Becker, J. Triimper: A&A 326, 682 (1997) W. Becker, J. Triimper: 
A&A 341 803 (1999) 

8. W. Bednarek, R.J. Protheroe: Phys.Rev.Lett. 79, 2616 (1997) 

9. M.C. Begelman, Z.-Y. Li: ApJ 426, 269 (1994) 

10. A.R. Bell: MNRAS 257, 493 (1992) 

A.R. Bell, S.G. Lucek: MNRAS 283, 1083 (1996) 

11. P. Bhattacharjee, G. Sigl: Phys.Rep. 327, 109 (2000) 

12. P.L. Biermann: Phys.Rev. D 51, 3450 (1995) 

13. O. Blanch, G. Blanchot, M. Bosman, et al.: 'The MAGIC telescope', Letter 
of Intent (1998) 

14. R.D. Blandford: preprint, [astro-ph/990602S; , (1999) 

15. P. Blasi, R.I. Epstein, A.V. Olinto: ApJ 533, L123 (2000) 

16. T. Bulik, B. Rudak, J. Dyks: MNRAS 317, 97 (2000) 

17. F. Camilo, et al.: ApJ 541, 367 (2000) 

18. K. Chen, M. Ruderman: ApJ 402, 264 (1993) 

19. K.S. Cheng, L. Zhang: Fundamentals of Cosmic Physics 20, 177 (1999) 
K.S. Cheng, L. Zhang: ApJ 515, 337 (1999) L. Zhang, K.S. Cheng: ApJ 
487, 370 (1997) 

20. K.S. Cheng, C. Ho, MA. Ruderman: ApJ 300, 500 (1986) K.S. Cheng, C. 
Ho, MA. Ruderman: ApJ 300, 522 (1986) 

21. F.V. Coroniti: ApJ 349, 538 (1990) 

22. J.W. Cronin: Rev. Mod. Phys. 71, S165 (1999) 

23. J.K. Daugherty, A.K. Harding: ApJ 252, 337 (1982) 

24. J.K. Daugherty, A.K. Harding: ApJ 273, 761 (1983) 

25. J.K. Daugherty, A.K. Harding: ApJ 309, 362(1986) 

26. J.K. Daugherty, A.K. Harding: ApJ 336, 861(1989) 

27. J.K. Daugherty, A.K. Harding: ApJ 429, 325 (1994) S.J. Sturner, CD. Der- 
mer: ApJ 420, L79 (1994) A.K. Harding, A.G. Muslimov: ApJ 500, 862 
(1998) 

28. CD. Dermer: ApJ 360, 197 (1990) 

29. J. Dyks, B. Rudak: A&A 362, 1004 (2000) 



30. J. Dyks, B Rudak, T. Bulik: Proc. of the 4th INTEGRAL Workshop, |astro- 
ph/001030l| , (2000) 

31. R.T. Emmering, RA. Chevalier: ApJ 321, 334 (1987) 

32. T. Erber: Rev.Mod.Phys. 38, 626 (1966) 

33. J.M. Fierro et al.: ApJ 447, 807 (1995) 

34. P. Goldreich, W.H. Julian: ApJ 157, 869 (1969) 

35. E.V. Gotthelf, G. Vasisht, M. Boylan-Kolchin, K. Torii: ApJ 542L, 37 (2000) 

36. E.M. de Gouveia Dal Pino, A. Lazarian: ApJ 535, L31 (2000) 

37. J.P. Halpern, F.Y.-H. Wang: ApJ 477, 905 (1997) 

38. A.K. Harding, M.G. Baring, P.L. Gonthier: ApJ 476, 246 (1997) M.G. Baring, 
A.K. Harding: to appear in ApJ 547, February 1 issue (2001) 



30 



Bronislaw Rudak 



39. A.K. Harding, I. Contopoulos, D. Kazanas: ApJ 525, L125 (1999) 

40. A.K. Harding, A.G. Muslimov: ApJ 508, 328 (1998) 

41. H. Herold: Phys.Rev.D 19, 2868 (1979) 

42. A.M. Hillas: A&A 22, 425 (1984) 

43. W.N. Johnson et al.: ApJS 86, 693 (1993) V. Schonfelder et al.: 
ApJS 86, 657 (1993) D.J. Thompson et al.: ApJS 86, 629 (1993) 

44. T. Kamae, T. Ohsugi, D.J. Thompson, K. Watanabe: Adv.Sp.Res. 25, 905 
(2000) 

45. C.F. Kennel, F.V. Coroniti: ApJ 283, 694 (1984) C.F. Kennel, F.V. Coroniti: 
ApJ 283, 710 (1984) 

46. L. Kuiper et al.: A&A 351, 119 (1999) 

47. L. Kuiper et al.: A&A 359, 615 (2000) 

48. S.G. Lucek, A.R. Bell: MNRAS 268, 581 (1994) 

49. F.E. Marshall, E.V. Gotthelf, W. Zhang, J. Middleditch, Q.D. Wang: 
ApJ 499L 179 (1998) 

50. S. Mereghetti: NATO ASI "The Neutron Star - Black Hole Connection" 
(1999) 

51. F.C. Michel: Theory of Neutron Star Magnetospheres. (The University of 
Chicago Press, Chicago 1991) P. Meszaros: High-Energy Radiation from Mag- 
netized Neutron Stars. (The University of Chicago Press, Chicago 1992) 

52. A.G. Muslimov, A.I. Tsygan: MNRAS 255, 61 (1992) 

53. H.I. Nel et al.: ApJ 465, 898 (1996) 

54. K. Nomoto, S. Tsuruta: ApJ 312, 711 (1987) J.M. Lattimer, K.A. Van Riper, 
M. Prakash, M. Prakash: ApJ 425, 802 (1994) 

55. H. Ogelman: 'X-Ray Observations of Cooling Neutron Stars'. In: NATO ASI, 
The Lives of Neutron Stars ed. M.A. Alpar, U. Kiziloglu, J. van Paradijs 
(Dordrecht: Kluwer 1995) pp. 101-120 

56. A.V. Olinto: preprint submitted to Elsevier Preprint, istro-ph/000200f: , 
(2000) 

57. R. Ong: ^stro-ph/0003014| (2000) 

58. G.G. Pavlov, Y.A. Shibanov, J. Ventura, V.E. Zavlin: A&A 289, 837(1994) 

59. A. Possenti, S. Mereghetti, M. Colpi: A&A 313, 565 (1996) 

60. M.J. Rees, J.E. Gunn: MNRAS 167, 1 (1974) 

61. B. Rudak, J. Dyks: MNRAS 295, 337 (1998) 

62. B. Rudak, J. Dyks: MNRAS 303, 477 (1999) 

63. Y. Saito et al.: ApJ 477 L37 (1997) Y. Saito, N. Kawai, T. Kamae, S. Shibata: 
'Search for X-ray pulsations from Rotation-Powered Pulsars with ASCA'. In: 
Neutron Stars and Pulsars: Thirty Years after the Discovery ed. N. Shibazaki, 
N. Kawai, S. Shibata et al. (Universal Academy Press: Tokyo 1998) pp. 295- 
298 

64. T. Sako, Y. Matsubara, Y. Muraki et al.: ApJ 537, 422 (2000) 

65. S. Shibata, N. Kawai, K. Tamura: 'Pulsar Nebulae as Calorimeters'. In: Neu- 
tron Stars and Pulsars: Thirty Years after the Discovery ed. N. Shibazaki, 
N. Kawai, S. Shibata et al. (Universal Academy Press: Tokyo 1998) pp. 457- 
464 

66. P. Slane: ApJ 437, 458 (1994) 

67. P. Slane, N. Lloyd: ApJ 452L, 115 (1995) 

68. V. Schonfelder et al.: A&AS 143, 145 (2000) 

69. PC. Schroeder, M.P. Ulmer, S.M. Matz et al.: ApJ 450, 784 (1995) 



Neutron Stars 31 



70. M.S. Strickman et al.: ApJ 460, 735 (1996) 

71. M.S. Strickman, A.K. Harding, O.C. de Jager: ApJ 524, 373 (1999) 

72. S.J. Sturner: ApJ 446, 292(1995) 

73. S.J. Sturner, CD. Dermer, F.C. Michel: ApJ 445 736 (1995) CD. Dermer, 
S.J. Sturner: ApJ 420, L75 (1994) 

74. P.A. Sturrock: ApJ 164, 529, (1971) 

75. M. Takahashi et al.: 'Detection of Pulsed X-ray Emission from The Fastest 
Millisecond Pulsar PSR B1937+21 with ASCA'. In: IAU Coll. 177, Pulsar 
Astronomy - 2000 and beyond ed. M. Kramer, N. Wex, R. Wielebinski (ASP 
Conference Series 202: 2000) pp.353-354 

76. J.H. Taylor, R N. Manchester, A.G. Lyne F. Camilo: 'The Princeton Pul- 
sar Database' (http://pulsar.princeton.edu) (1995) J.H. Taylor, R.N. Manch- 
ester, A.G. Lyne: ApJS 88, 529 (1993) 

77. D.J. Thompson, A.K. Harding, W. Hermsen, M.P. Ulmer: 'Gamma-Ray Pul- 
sars: the Compton Observatory Contribution to the Study of Isolated Neutron 
Stars'. In: Fourth Compton Symposium ed. by CD. Dermer, M.S. Strickman, 
J.D. Kurfess (New York: AIP 1997) 410 

78. D.J. Thompson et al.: ApJ 516, 297, (1999) 

79. A. Venkatesan, M.C Miller, A.V. Olinto: ApJ 484, 323 (1997) 

80. H. Volk, HEGRA Collaboration: 'HEGRA Observations of Galactic Sources'. 
In: GeV-TeV Gamma Ray Astrophysics Workshop : towards a major atmo- 
spheric Cherenkov detector VI ed. B.L. Dingus, M.H. Salamon, D.B. Kieda: 
AIP Conference Proceedings, Vol. 515., p. 197 (New York: AIP 2000) F. Aha- 
ronian et al.: A&A 346, 913 (1999) R.W. Lessard et al.: ApJ 531L, 942 
(2000) 

81. T.C Weekes: 'VHE Astronomy Before the New Millennium'. In: GeV- 
TeV Gamma Ray Astrophysics Workshop : towards a major atmospheric 
Cherenkov detector VI ed. B.L. Dingus, M.H. Salamon, D.B. Kieda: AIP 
Conference Proceedings, Vol. 515., p.3 (New York: AIP 2000) M. Catanese, 
T.C. Weekes: PASP 111, 1193 (1999) 

82. T.C. Weekes, S.M. Bradbury, I.H. Bond et al.: 'The Very Energetic Radiation 
Imaging Telescope Array System (VERITAS)', preprint (2000) 

83. XY. Xia, G.J. Qiao, X.J. Wu, Y.Q. Hou: A&A 152, 93 (1985) H.K. Chang: 
A&A 301, 456 (1995) 

84. M.D. Young, R.N. Manchester, S. Johnston: Nature 400, 848 (1999) 

85. B. Zhang, A.K. Harding: ApJ 532, 1150 (2000) 

86. B. Zhang, A.K. Harding: ApJ 535, L51 (2000) 



